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One-dimensional non-equilibrium kinetic Ising models evolving under the competing effect of 
spin flips at zero temperature and nearest neighbour spin exchanges exhibiting a parity-conserving 
(PC) phase transition on the level of kinks are now further investigated, numerically, from the point 
of view of the underlying spin system. Critical exponents characterising its statics and dynamics are 
reported. It is found that the influence of the PC transition on the critical exponents of the spins is 
■ strong and the origin of drastic changes as compared to the Glauber-Ising case can be traced back 

to the hyperscaling law stemming from directed percolation. Effect of an external magnetic field, 
leading to directed percolation type behaviour on the level of kinks, is also studied, mainly via the 
generalized mean field approximation. 
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I. INTRODUCTION 



Nonequilibrium phase transitions have attracted great interest lately. A variety of sytems studied seem to belong 
to the universality class of directed percolation (DP) K-m The DP universality class is very robust, among others 
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the transitions in branching annihilating random walk (BARW) models with an odd number of offsprings in the 
process A — > A + nA belong to it || . Numerical studies by Grassberger et al. HQ of probabilistic cellular automata 
^\ , models in one dimension involving the processes k — > 3fc and 2k — > (k stands for kink) have revealed , however, a 
■ new universality class of dynamic phase transitions. Both time-dependent and steady-state simulations have resulted 
non-DP values for the relevant critical exponents. This so called parity-conserving (PC) phase transition has been 
found since in a variety of models. The n — 4 BARW model has been studied in the greatest detail and accuracy 
by Jensen . In a previous paper [^j one of the the authors has introduced a family of nonequilibrium kinetic Ising 
models (NEKIM) showing the same phenomenon, while quite recently, with an appropriate modification of the original 
BARW model (suggested already in ref. H), the n — 2 BARW model has been shown to exhibit the PC transition 
as well JTo[ . The two-component interacting monomer-dimer model introduced by Kim and Park ]Tl|-p^] represents 
a more complex sytem with a PC-type phase transition. 

Formerly, combinations of Glauber juj and Kawasaki [ Jl 5 | kinetics were introduced with the aim of investigating 
temperature-driven nonequilibrium phase transitions |l6| , |l7|] . In rcf. |J , however, spin-flip kinetics was taken at T = 
alternatingly with random nearest-neighbour spin-exchanges ( Kawasaki-exchange at T — oo). The mean- field limit 
of this model together with results of the generalised mean-field theory (GMF) have been presented in |Is|] . 

In the present paper we further investigate NEKIM but now from the point of view of the underlying ID spin system 
with the aim of getting some more insight into the nature of the PC transition. Results of computer simulations are 
presented for different critical exponents connected with the structure factor. The new computational technique 
based on short-time dynamics and finite size scaling, as introduced and applied by Li,Schiilke and Zheng |lj],|2(| for 
calculating moments of the magnetisation and the time dependent Binder cumulant, is also applied here besides the 
usual finite size scaling (FSS) and time dependent simulations. In all our numerical studies the initial state is random 
with zero magnetisation and concerning finite size effects the use of antiperiodic boundary conditions instead of the 
usual periodic ones will prove to be essential. 

In the field of domain- growth kinetics it has long been accepted that the scaling exponent of L(t), the characteristic 
domain size, is equal to 1/2 if the order parameter is non-conserving. Now let us restrict ourselves to a ID Ising spin 
chain of length L and define the structure factor as usual : S(0, t) = L[< M 2 > — < M > 2 ], M = j- ^ Sj, (s; = ±1). 
If the conditions of validity of scaling are fulfilled |H| then 

5(0, i)) oc [L(t)] d ; L{t)cxt x (1) 

where d = 1 now and x = 1/2. An other usually considered quantity is the excess energy of a monodomain sample at 
the temperature of quench, which in our case is proportional to the kink density n(t) = < J2i 5(1 — s i s i+i) >• 
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with y = 1/2 in the Glauber-Ising case, expressing the well-known dependence on time of annihilating random walk. 

It is also well-known, that the parameter pr — e ~ 4J / kT can be regarded as the quantity measuring the deviation 
from the critical temperature T — of the ID Ising model, and the usual (equilibrium) critical exponents can be 
defined as powers of pr- From exact solution the critical exponents of the spin- susceptibility, coherence length and 
magnetisation are known to be 7 = v = 1/2, (3 = 0, respectively. Fisher's static scaling law 7 = dv — 2(3 is valid. 
Moreover, the dynamic critical exponent Z is equal to 2. 

Instead of Glauber kinetics, let us apply now the nonequilibrium NEKIM kinetics, and investigate the Ising system 
at and in the immediate neighbourhood of one of its PC points. On the level of kinks the PC point, separating the 
active and absorbing phases, is a second order transition point , while from the point of view of spins the absorbing 
phase consists of a multitude of ID Ising critical points, which ends at the PC point. We are interested in the 
behaviour of the spin system at this end-point. 

The result is the following: the critical kink-dynamics has strong influence on the spin-kinetics and even on its 
statics. Domain growth is governed by criticality: x=l/Z, but the dynamic exponent Z changes from Z = 2 to 
Z = 1.75. As to statics, 7 = v = .AAA while (3 = and Fisher's scaling law remain valid. 

More spectacular is the change of the above y-exponent: y = x/2;x = .57 and similarly the exponent a n of 
the finite-temperature kink-density: limt-yoo n(t, pr) oc pr an , which is 1/2 in the Glauber-Ising limit, decreases to 
a n = 7/2 = .222. We shall argue that this factor of 2 between magnetic- and kink-exponents has its origin in the 
hyperscaling law, introduced for DP by Grassberger and de la Torre |1| , which connects exponents of time-dependent 
kink density and cluster-size. 

These numerical results will be reported in detail in the following. Moreover, we shall present further results, mainly 
from the side of applying GMF approximations, by introducing a magnetic field-term into the spin- flip probability, 
which causes the PC transition to become DP-like [ p^| . 



II. THE MODEL 



The model we will investigate here is a one-dimensional kinetic Ising model evolving by a combined spin-flip and 
spin-exchange dynamics as described in |1. The spin-flip transition rate in one-dimension for spin s,; sitting at site i 
is: 

r / 7 \ 

Wi = -(1 + 5si-is i+ i) (1 - -s i (s l ^ 1 + s i+1 )j , (3) 

where 7 = tanh 2 J/fcT (J denoting the coupling constant in the Ising Hamiltonian) , T and 6 are further parameters. 
While ia H T = (7 = 1) has been taken, we shall consider now finite temperature effects, too. Instead of 7 the 
parameter px will be used in the following. The three independent rates: 

r 1 pt 

Windif = ttC 1 ~ <*)) w op P o = r(i + 8)— — , w same = r(i + 6)— — , (4) 

2 1 + pr I+Pt 

where the subscripts of w refer to the three possible neighbourhoods of a given spin, are responsible - on the level of 
domain walls - for random walk, annihilation and pairwise creation (inside of a domain) of kinks, respectively. 
The other ingredient of NEKIM has been a spin-exchange transition rate of neighbouring spins ( the Kawasaki fill 
rate at T = 00): 

Wu+l = ~ SiS i+1 ], (5) 

where p ex is the probability of spin exchange. 

Spin-flip and spin-exchange have been applied alternatingly at each time step, the spin-flip part has been applied 
using two-sublattice updating, while making L MC attempts at random ( L denotes the size of the chain) has been 
counted as one time-step of exchange updating. In this system, at T = 0, PC type phase transition takes place. In || 
we have started from a random initial state and spotted the phase boundary in the {5, pex) plane. In the following 
we will choose a typical point on this phase diagram and make simulations at and around this point, fixing T and p ex 
and changing only 8. The parameters chosen are: V = .35, pex = .3, 5 C = —.395(2). We note here that T appeared as 
1 /r in ref . |£j] and 5 C has a more accurate value now than previously. 
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III. SCALING FORMS AND LAWS 



In the following we will be interested in two quantities characterising the behaviour of the NEKIM, namely the 
structure factor or spin-susceptibility and the kink-density under the conditions of a quench. Thus, in contrast to 
the usually considered evolution from a pair of kinks, we will restrict ourselves to completely random initial states ( 
T = oo, M (0) = 0) and follow the development of the system via the rules described in the previous section. With pr 
defined above and e =\ S — S c , where S < is the parameter which drives the phase transition in the present case, 
the scaling forms for S and n are as follows: 

t~h ti £ £i 

S(p T> e,L,t)=t*f(— (6) 



n(p T ,e,L,t) = t-yg(— ,— ,|,^) (7) 

with £ = pt~ v , £_L = £~ u± ■ The exponents connected with e have been written using the notation of directed 
percolation, while those related to the temperature factor px are written in the notation of equilibrium Ising system. 
Z and Z c are the respective dynamic critical exponents and we allow for the possibility that they differ (though this 
will turn out not to be the case). 



A. The large L case 

Let us hrst take the limit L — » oo. Then the dependences on £/L and £j_/L can be neglected in eqs(||),(0). If 
furthermore £ — > oo, £^ — > oo the forms valid at the critical point e = 0, pr = are obtained: S c (t) oc t x and 
n c {t) oc t~ y . Fig.l. shows S c (t) oc t x , with the result x — .570(1). The same statistics has led now for the exponent 
of the kink density to the value y — .285(1) (denoted a in [9]). 




FIG. 1. The sructure factor as a function of time at a typical phase transition point of the NEKIM phase diagram. L = 128 
and the number of independent initial configurations with zero magnetisation was 10 s . Here tm ax = 200 and finite size effects 
start to set in at about t = 1000. The straight line is fit with x = .570. 
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We have also made simulations for the same quantities with L = 8000 up to t — 3 x 10 4 and averaging over 2000 
initial states, with similar result but with higher error and more than 100 times bigger computer time. It is worth 
noticing, that within error x — 2y. 

Keeping px — 0, near S c in the active phase eq.(Q) gives, for t — > oo the order parameter n(e) cx e' 3 " with the scaling 
law § 

p n = v^Zy. (8) 

A similar scaling relation can be obtained from (^), which yields for t — » oo : S(e) oc e -6 with 

6 = u ± Zx. (9) 

The divergence of the spin-susceptibility as a function of e is understandable: in the subcritical regime (i.e. for 
S \<\ S c | ) it is infinite ( T — > 0, t — > oo), because the whole subcritical region is a plane of Id critical (Ising) 
points. In order to get directly we have made simulations for S(e,t) around the chosen PC point in the interval 
e = .02 — .13 as shown on Fig. 2. 
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t 



FIG. 2. The structure factor S(e, t) in the active phase in the vicinity of the PC transition point for the following values 
of e: .13, .11, .09, .075, .06, .04, .03 (from bottom to top). L was varied between 2000 and 8000 to avoid finite size effects 
before saturation sets in. Number of independent states in averaging: typically 1000. 

The (time-averaged) saturation values of S yield the exponent = 1.9(1) . From the same runs as above we have 
obtained now (3 n = .88(4). Thus within error O = 2/3„, which relation also follows from eqs. (||) and (||) with x = 2y. 
It is worth noting here that on the basis of the divergence of the spin susceptibility, S(e), the PC transition point (as 
endpoint of a line of first order transitions) can be found, without any reference to kinks. 

Taking now e = and keeping pr finite in the limit t — > oo we get from eq. (0) n(pT) oc pt u " and the scaling 
relation 

a n = yvZ c . (10) 

Fig. 3 shows n(t,pr) as a function of t at different values of pr at e = 0; the level-off values could be fitted with 
a n = .222(5). This result is in accord with the value for the n = 4 BARW reported by Jensen S (1/Sh in his 
notation). 

Similarly, taking eq.(||) at e = and in the limit t — > oo the spin-susceptibility arises: x 00 Pt ~ 7 together with the 
scaling relation 

7 = xvZ c . (11) 
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On Fig. 4 S(px, e = 0,t) is plotted for different values of pr- 
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FIG. 3. The kink-density as a function of temperature and time at the PC point. The curves are for 
pr = 1CT 1 , 5.0 x 1(T 2 , 1(T 2 , 5.0 x 10 -3 , 10 -3 , 5.0 x 10~ 4 from top to bottom. L = 2000 and the number of samples 
in averaging typically 2 x 10 4 . 




FIG. 4. The structure factor S(pr,t) for different values of pr- Details of simulation as of Fig. 3 except that pr is decreasing 
from bottom to top here. 

The temperature-dependent simulations have been performed in the range px = 1.10 -2 — 5.10 -5 . The level-off 
values seen on Fig. 5, where the data for the two highest temperatures ( pr = 1.10 -1 , 5.10 -2 ) were discarded , yield 
the exponent of the spin-susceptibility at the PC point as 7 = .445(5). This value is in accord, within error, with the 



scaling laws (10) and (O) which predict , using the relation x = 2y found above, 7 = 2a n = -444 . 
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FIG. 5. The saturation values of Fig. 4 time-averaged, on a double-logarithmic scale. Data for pr = 10 4 , 5 x 10 5 , not 
shown on Fig. 4, have been included. Straight line: fit of data points with 7 = .445 

It is to be noted, that due to the non-self-averaging property of the structure factor |^2| , all exponents connected with 
this quantity have much larger statistical errors than the ones connected with the kink density, and the same applies 
to the time-dependent simulations. This explains the large fluctuations exhibited on Fig. 2 and Fig. 4 in comparison 
with e.g. Fig. 3. 



B. Static finite size scaling 

Finite size scaling will be used to find numerical values for some more exponents. The static FSS limit requires to 
take first the limit t — > 00 and to suppose L <c £, Let us consider now the two possible orders of limits to reach 
the PC point: a), pt — 0, e — > and b). e = 0, px — ► 0. 

Case a). 

Eqs (H) and (0) lead to the expressions: ]im e -^oS(e, L) = L~ f"(-^-) cx L~ and 

lim^Q n(e , L) = L v± 9 oc L "x , respectively. 
Case b). 

From eqs (^|) and (Q) we get now lim PT ^ S(pr, L) = L~ f"(j) cx L~ and 
lim PT ^o n (pr, L) = L~~i?~g"(j) oc L^~^, respectively. 

Within error we have not found any (numerical) evidence against the supposition that the order of limits e — > 
and pr — > were interchangeable, the same fixed point seems to be reached in both cases. ( Actually this fact has 
expressed itself already in our finding, according to which the time-exponents x and y are also the same independently 
from where we approach the limit e = 0, px = ; the difference should have been seen in the preasymptotic time- 
dependence. E.g. for very small values of pt at e = x = .57, y = .28 was always clearly seen in the early-time 
behaviour of S(t) and n(t), resp.). 

Thus the above relations lead to the following scaling equalities: 

7 = Pn _ On 

v u± 1 l/± V 

In numerical simulations with random initial states usually periodic boundary conditions (pbc) are supposed. The 
fact that pbc allow only an even number of kinks does not matter except under the conditions of FSS close to the 
PC point as eventually for all samples the ordering becomes perfect (depletion of kinks) and trivially J = = 1. 
Because of the same reason it is not possible to find the L-dependence of the kink density, which breaks down as a 
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function of time. The proper procedure is to use antiperiodic boundary conditions (apbc), which choice allows only 
an odd number of kinks (i.e. all samples are surviving) provided the updating procedure is carefully done. As in the 
course of dynamic FSS one reaches fixed points, which coincide with the limiting values for t — > oo of S and n above, 
we will cite our results in the next section. 

As we shall see 2- = 1 is valid also in case of antiperiodic boundary conditions, which value together with eq.(|Tl|) 
gives Z c = 1/x. Moreover, we shall arrive at = 1/2, which, using eqs. (8) and (9), gives ^- = 1 and leads 
to Z = 1/x. Thus Z c = Z, the two critical dynamical exponents coincide, as anticipated. It is worth noting that 
Z = type relation between dynamic critical exponent and domain growth exponent has been found earlier for 
various lei equilibrium Ising systems j2^| . 



C. Early-time dynamic Monte-Carlo method 

In case of systems quenched to their critical temperature p4[ universality and scaling may appear in a quite early 
stage of time evolution, far from equilibrium, where £ is still small. Based on the scaling relation for such early time 
intervals, a new way for measuring static and dynamic exponents has been proposed p9| , |20| |. Now we apply this 
method to get critical exponents for a nonequilibrium phase transition. 

Following we shall suppose the following relation to hold for the k-th moment of the magnetisation near the 

critical point of the ID spin system: 

M (k \t, PT ,L) = b~ h ^M {k \b- z t,b 1/u PT ,b- 1 L) (13) 

where zero initial magnetisation has been considered and b is a rescaling factor (b = 2 will be chosen). After generating 
randomly an initial configuration, the system is let to evolve according to the nonequilibrium kinetic rule of Section 
II at pr — 0. ( We could have included e in eq.(|l3|) as well, but we will restrict ourselves to the case e — 0, pr = 0, 
so it is of no importance here). Average is taken over the initial configurations with zero magnetisation. In order to 
get sufficiently good statistics averaging has to be performed over very many (10 5 — 10 6 ) independent initial states, 
as emphasised by Li et al. Jl5^(| who applied the method to the 2d Ising model. To get the dynamical exponent Z 
with great accuracy, they proposed calculation of the time-dependent Binder cumulant: 

m '^ L) = 1 -^f (14) 

which behaves at pr = as 

U(t, 0, L) = U{b- Z t, 0, L/b) (15) 
Similarly to eq. (|l3|), n(t, L) can be supposed to behave according to the relation: 

n{t, L) = b~^Tn(b- z t, b~ l L) (16) 

Besides the triviality of the fixed points reached (M^* = 1, U* = 2/3) the pbc case still leads - by proper fitting 
- to the value of Z and (3. In case of M< 2 \ using eq. (ft3|), collapsing of curves for different values of L has lead to 
/3 = 0.00(2), Z — 1.75(2); see Fig. 6. In case of U (t, 0, L), which contains only Z as fitting parameter and thus provides 
it directly, collapsing of curves could be achieved again with Z = 1.75(1). Typically 10 5 — 3.10 5 averages have been 
performed. 
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FIG. 6. Collapsing of curves for M 1 - 2 ' with periodic boundary conditions on a double-logarithmic scale. A -s are rescaled 
data while • -s are original ones. From top to bottom: L=16(A)-32(»), 32(A)-64(«), 64(A)-128(»), 128(A)-256(»). Number 
of samples with random initial states in the averaging was typically 10 5 



With antiperiodic boundary conditions collapsing of curves has led us to the same results as above concerning the 
exponents f3 and Z . Under these conditions, however, important additional information arises from best fit to time 
averaged saturation values. Concerning n{t,L), eq.(|l^), the collapse of curves is illustrated with L = 100,200 on 
Fig.7. 




t 

FIG. 7. Collapsing of curves for n(t,L) according to eq.(16). A-s: L=100, rescaled; »-s: L=200. Number of samples in the 
averaging: 10 5 . Boundary condition: antiperiodic, scale: double- logarithmic. 



The curves scale together with Z = 1.75(1) and £j = .52(2). On the other hand, fitting the time-averaged 
saturation values for L — 50, 64, 100, 128, 200, 256 gives for the kink-density n(L) = .66 x L~ AS( - 2 l Thus we conclude 



8 



|a.(= 2a.) = .50(2). This is in accord with the result of Jensen for the n = 4 BARW ||. 
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FIG. 8. M^ 2 \t,L) with apbc. Rescaling with & = 2, /3 = and Z = 1.75 has resulted in data points marked by A -s. 
Collapsing of curves is shown -from top to bottom- for L=32(A)- 64(»), 64(A)-128(«), 128(A)— 256(»). Number of samples 
in averaging: 10 5 — 3 x 10 5 




100 1000 
t 

FIG. 9. U(t,L) with apbc for L=16(A)-32(«), 32(A)-64(»), 64(A)-128(«) (from top to bottom). Best collapse of curves 
again wit Z = 1.75. Number of samples in averaging : as for Fig.8. 

On Figs. 8 and 9., respectively, M^> and U are seen with antiperiodic boundary conditions for some values of L. 
It is apparent, that U is much more sensitive: here even much more samples in averaging would have been needed to 
smooth out the curves. The time-averaged saturation values for L = 50, 64, 100, 128, 200, 256 lead to S{L) = .26£' 99 W 
(Fig.10), thus £(= = .99(1). 
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FIG. 10. Time-averaged saturation values S(L) of M^ 2 \t,L) shown on a double-logarithmic scale. Straight line: power-law 
fit with j/v = .99. Averaging has been performed over 10 5 — 3 x 10 5 independent random initial states with apbc. 

The difference as compared with the pbc case shows up in the prefactor ( .26(1) instead of 1.0 ). For the Binder 
cumulant U* = .32(1), again differing from the pbc value of 2/3. For the sake of comparison, we have carried out 

the same kind of simulations for the Ising-Glauber case (p ex = 0, S — 0) with apbc, and obtained M^* — .333(3), 
n* gi = 1.00(5)/i, U*gi = -40(1) . Here Z = 2 and (3 = have given the best fit and similar result can be expected 
in the whole absorbing phase. It is worth mentioning here, that these fixed-point values can be derived exactly for 
the Glauber case with the result: M^ } * = 1/3, n* G i = 1/L, U* Gi = .4 §|. 



IV. SYMMETRY-BREAKING FIELD 



It is by now well established that the PC transition has non-DP critical exponents, because of the modulo 2 
conservation law. Park and Park Jl3| have introduced symmetry breaking external field in case of the interacting 
monomer-dimer model and showed that the DP universality class is recovered if one of the absorbing states is singled 
out. They also have mentioned of having similar, though - to our knowledge - not yet published data for NEKIM 
and Grassberger's automata. We have investigated the effect of an external magnetic field H on the NEKIM model 
with simulations and with the help of the Generalized Mean-Field (GMF) technique as well, both confirming the DP 
behaviour. The transition probabilities of NEKIM, as given in Sect. II., are modified in the presence of en external 
magnetic field H as: 

Windif = W indl f {I -hSi), (17) 

oppo (1 - hSi), (18) 
h = th(^). (19) 



w oppo ' ' ' ' 



Here we shall restrict ourselves to pt = 0. Fig. 11 shows the phase diagram of NEKIM in the (h,S) plane, starting at 
the reference PC point for h — used in this paper (S c = —.395). 
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FIG. 11. Phase diagram of NEKIM in the (h, S) plane in the presence of an external magnetic field. Parameters of the 
transition probabilities: V — .35,p ex = -3, see Section II. Naturally, the phase diagram can be drawn symmetrically for negative 
values of ft, as well. 

We have applied only random initial state simulations to find the points of the line of phase transitions (critical 
exponents: y=.17(2), (3 n — .26(2)). It is seen, that with increasing field strength the critical point is shifted to more 
and more negative values of <5. 

As to the treatment of the model with the use of the GMF technique, the details of which have been explained 
in [2fj-p£|, in the field-free case we have already obtained estimates for the critical point and the effect of particle 
exchange earlier |L8| . By applying the Coherent Anomaly extrapolation method (CAM) J2^] one can extrapolate the 
critical exponents of the true singular behavior and in [ p"8[ we could give a rough estimate for the exponent f3 n of the 
PC transition based on N < 6-th order cluster GMF calculation. Now we extend the method for the determination 
of the exponent of the order-parameter fluctuation as well : 

Xn {e) = L(<n 2 > - <n > 2 ) ~ e^* . (20) 

The GMF equations have been set up for the steady states of NEKIM in the presence of the H-field. The A^-block 
probabilities were determined as the numerical solution of the GMF equations for N — 1, 6. The traditional 
mean- field solution (N = 1) results in stable solutions for the magnetization : 

M = -j, if 5 < and h 2 /5 2 < 1 (21) 

M = sgn(h), otherwise. (22) 

and for the kink-concentration : 

ra=±(l-(§) 2 ), if 5<0 and h 2 /S 2 < 1 (23) 

n = 0, otherwise. (24) 

For TV > 1 the solutions can be found numerically only. By increasing the order of approximation the critical 
point estimates 5 C (N) shift to more negative values similarly to the H = case. The limN^oo S c (h) values have been 
determined with quadratic extrapolation in case of h = 0.01,0.05,0.08,0.1. The resulting curves for n(S) and Xn(S) 
are shown on Figures 12. and 13., respectively for the case of h = .1 in different orders N of the GMF approximation. 
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FIG. 12. The kink density in the neighbourhood of the critical point 5 c (h) for h — .1. The curves from right to left 



correspond to TV = 1, 6 (level of GMF calculation) 
to be able to extract CAM anomaly coefficients. 



The points have been determined with resolution of 10 in 6, in order 
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FIG. 13. The same as Fig.12 but for the second moment of the kink density. It is apparent also here that the value of 8 c (h), 
being equal to .70(1) for h = .1 according to simulations, is fairly well approximated by the GMF-results for N — 5,6. 



Naturally, these curves exhibit mean-field type singularity at the critical point 

n(N) ~-p n (8/6 c (N)-l)f , >" 
Xn(N)~x n (N)(5/6 c (N)-l)-™*, 



(25) 
(26) 



with Pmf = 1 and -fMF = — 1- According to the CAM (based on scaling) the critical exponents of the true singular 
behavior can be obtained via the scaling behavior of anomaly factors: 



n(N) 



(27) 
(28) 



where we have used the A = (S C /6 C (N) — 5 C (N)/8 C ) invariant variable instead of e, that was introduced to make the 
CAM results independent of using S c or 1/8 C coupling ( fl30f ). Since the level of the GMF calculation, what we could 
solve is N < 6, we have taken into account correction to scaling, and determined the true exponents with non-linear 
fitting form : 



n ( N) = a A /3 '*~ /3mf + b A^-Pmf+i 



(29) 
(30) 
(31) 
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where a and b are coefficients to be varied. The results for various h-s are shown in table below. For h = 0.0 we could 
not determine the 7„ exponent, because the low level GMF calculations resulted in discontinuous phase transition 
solutions - what we can not use in the CAM extrapolation - and so we had too few data points to achieve a stable 
non-linear fitting. Higher order GMF solutions would help, but that requires the solution of a non-linear set of 
equations with more than 72 independent variables. This problem does not occur for h ^ 0; the above results - being 
based on all N = 1..6 point approximations - arc fairly stable. 



TABLE I. CAM calculation results 



h 


0.0 


0.01 


0.05 


0.08 


0.1 


DP 




1.0 


0.281 


0.270 


0.258 


0.285 


0.2767(4) 


7" 




0.674 


0.428 


0.622 


0.551 


0.5438(13) 
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V. SUMMARY AND CONCLUSIONS 



Time dependent simulations, FSS and dynamic early-time MC method have been applied here to investigate the 
behaviour at and in the vicinity of a PC transition of the nonequilibrium kinetic Ising model with the aim to complete 
earlier results. Emphasis has been put now on the critical properties of the Id spin system underlying kinks. In 
this way we have arrived at a more-or-less complete picture of the effect the PC transition exerts on the statics and 
dynamics of the ID Ising spin phase transition: we have found Z — 1.75 (instead of Z = 2.0); 7 = v = .444 (instead 
of 1/2). The PC point is endpoint of a line of first order phase transitions ( by keeping p ex and T fixed and changing 
8 through negative values to S c ); , where still (3 — holds and also Fisher's scaling law: 7 = dv — 2/3 is valid. The PC 
point has been approached from two directions: from the active phase by changing e = 8 C — 8 and from the direction of 
finite temperatures by varying pt — e~ 4J / kT . The second moment of the magnetization ( structure factor ) provides 
the 'magnetic' exponents x, 6,7, which have been found - within error - to be twice as big as the corresponding kink 
exponents, y,{3 n , a n . The cause of this factor of two must lie in the nature of the active phase, of course, and it is 
sufficient to understand x — 2y (time-dependent exponents at PC) as the rest follows from scaling relations. 

The idea is to recognize that there are two characteristic growth lengths at PC which have to have the same critical 
exponents. Namely, the magnetic one L(t) oc t x (see Introduction) and the cluster size defined through the square-root 
of < R 2 (t) >oc t z . The latter one is obtained by starting either from two neighbouring kink initial states (see e.g. 

), or from a single kink jTjIj, while the magnetic domains grow in the quenching situation i.e. from random initial 
states (T = 00 — > T = 0). Both length-exponents are, however connected with Z, the dynamic critical exponent, 
since at PC the only dominant length is the ( time- dependent) correlation length, § = ^ has been shown to follow 
from scaling in [jf) and ]3l| ] for one-kink and two-kink initial states, respectively. Presently we have found x = l/Z, 
thus x — z/2 follows. 

Exponents of kink density and cluster growth are connected by a hyperscaling relation first established by Grass- 
berger and de la Torre jl) for the directed percolation transition. In the same form it does not apply to PC transition 
|^|,[52j, where dependence on the initial state (one or two kinks) manifests itself in two cluster-growth quantities: the 
kink- number N(t) oc t v and the survival probability P(t) oc t~~ s . (This 8 has, of course, nothing to do with the 
parameter 8 of NEKIM) In case of odd number of kinks (when 8 — 0, as all samples survive) hyperscaling reads [Q: 
y + rj = z/2 As both y and 77 are exponents of time dependent kink density at PC, they should agree numerically 
if scaling is valid with a single (time-dependent) characteristic length, and this is indeed the case (see e.g. || for 
simulation data). Consequently 2y = z/2 — x follows from the hyperscaling relation. 

For an even number of kinks (two-kink situation) the relation with the random initial-state situation investigated 
here can be established as follows. The local kink density is defined as p(x, t) oc t v ~ z / 2 f(x 2 /t z ) which divided by P(t) 
leads to our n(t) of eq.(|^) at PC: n(t) oc t~ y with y = —rj — 8 + z/2. Using the hyperscaling relation valid in the 
two-kink case, namely 28 + r\ = z/2, y = 8 follows. Thus the hyperscaling relation can be written as 2y + r\ — z/2. 
Taking into account the result first obtained numerically by Jensen ]8| that i] = for the PC transition when starting 
with two kinks, 2y = z/2 — x holds, again. 

In this way the factor of two between magnetic and kink exponents found in this paper could be explained as 
following from scaling. 

In the last chapter we have introduced a magnetic field into the NEKIM transition rates and investigated its effect 
mainly in the framework of generalized mean field approximation. By going up to N=6 order cluster approximation, 
the expectation that the universality class of the phase transition turns into DP-type Jl3| has been given support: 
we have found values for the exponents of the kink density and its second moment which are very close to the corre- 
sponding DP values. 
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